Fractional ^-scaling for quantum kicked rotors without cantori 
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Previous studies of quantum delta-kicked rotors have found momentum probability distributions 
with a typical width (localization length L) characterized by fractional 7i-scaling, ie L h?^^ 
in regimes and phase-space regions close to 'golden-ratio' cantori. In contrast, in typical chaotic 
regimes, the scaling is integer, L ^ . Here we consider a generic variant of the kicked rotor, 
the random-pair-kicked particle (RP-KP), obtained by randomizing the phases every second kick; it 
has no KAM mixed phase-space structures, like golden-ratio cantori, at all. Our unexpected finding 
is that, over comparable phase-space regions, it also has fractional scaling, but L ~ h~^^^. A 
semiclassical analysis indicates that the h^^^ scaling here is of quantum origin and is not a signature 
of classical cantori. 

PACS numbers: 05.45Mt,32.80.Pj 



A notable recent achievement of atom optics has been 
the realization of a well-known paradigm of quantum 
chaos, the quantum '(5-kicked particle' (DKP) [l|, 
Laser-cooled atoms in a pulsed periodic optical potential 
can be effectively modeled theoretically by the Hamilto- 
nian H = p"^ /2 — kcosxJ2n 5{t — nT) provided the pulses 
are sufficiently short. This system has been extensively 
investigated in numerous theoretical (see eg [1, 01) and 
experimental [5] studies by several groups worldwide. In 
the chaotic regime {kT 3> 1) the momentum distribu- 
tions N{p, t) evolve into a final, time-independent distri- 
bution of exponential form: N{p, t — > oo) ~ exp(— 
L, the momentum localization length, has well-known in- 
teger scaling properties L ^ K^Ti^^ . This effect, termed 
'dynamical localization (DL)', has been well-studied both 
experimentally and theoretically. 

But an important exception, with L ~ h?^^ , was found 
in the seminal theoretical study of quantum behavior 
in the vicinity of so-called 'golden-cantori' in 0. At 
a critical value of kT w 0.97, the last classical barrier 
(KAM torus) which impedes chaotic diffusion is broken. 
What remains are fractal partial barriers, termed can- 
tori, situated at momenta p ~ 2-kR and p « 2tt{R — 1), 
(or integer multiples thereof) where R is the golden ra- 
tio. A subsequent study [7| suggested that a positive 
exponent L ~ Tz."*""^ was associated with tunnelling trans- 
port (favoured by increasing K) while a negative expo- 
nent L ~ was associated with dynamical localization 
where transport increases as S — > 0. In Q it was found 
that the sign of the scaling exponent can change from 
negative to positive as the dominant transport mecha- 
nism changes from tunnelling to dynamical localization. 
A study of the classical phase-space scaling near these 
golden tori Q identified two characteristic scaling ex- 



ponents, a « 0.65 and a « 0.75. Since the fractional 
L ^ fi' ^ h^^^ quantum scaling was found in regions 
close to p ~ 27ri?, this behavior was attributed in 6] to 
the smaller of the classical exponents. 

A recent study [l^ of a closely related system, the 
double- kicked particle (2-DKP), found fractional scaling 
of momentum distributions with L ~ ?i~'^'^^,ie charac- 
terized by one of the classical golden ratio exponents - 
and in phase-space regions corresponding quite closely 
to those of [6|]- it was argued that this too was evidence 
of the quantum signature of the golden cantori. As the 2- 
DKP has already been experimentally realized with cold 
atoms Q and its fractional scaling occurs over a much 
broader range of p than for the usual DKP, this feature 
of the atomic momentum distributions is very amenable 
to experimental study. 

However, here, we introduce a system which is ran- 
dom, but closely related to both the DKP and 2-DKP: 
the randomized-pair DKP or RP-DKP. It is obtained 
as a limit as the phases between consecutive kick-pairs 
become completely random. This system has no mixed 
phase-space regions or KAM cantori at all. This study, 
over a range of smaller Ti values than [13] found the ran- 
dom system has exactly the same scaling properties as 
the mixed-phase-space 2-DKP. Both have L ~ h~^^^ for 
small h; in a slightly larger exponent was found be- 
cause some non-asymptotic values of Ti were included. 
This unexpected result is the key finding of our paper. A 
semiclassical analysis shows that the scaling with Tc"^^^ 
arises naturally from the quantum unitary matrix and is 
not a signature of golden cantori. We discuss possible 
implications for the ordinary DKP. 

The Hamiltonian of the 2-DKP is given by = 
^-kcosxY^^W -]<iT)^E(t' -NST^t)] where e < 



T is a small time interval. In effect, the particles are 
exposed to a sequence of pairs of closely spaced J-kicks. 
One could now adopt the usual procedure and choose T 
to define the time unit and thence derive a classical map 
with a stochasticity parameter K = kT and momentum 
units p = P'T. However, we show below that, for the 2- 
DKP, it is the small timescale e, rather than the period T, 
which provides the natural unit of time. Hence, rescaling 



time as t'/e 
map: 



t, and P'e p, we obtain our 2-DKP 



PN+i =pn + sinxjv; Pn+2 = Pn+i + sinxjv+i 

XN+l = XN + PN+l] XN+2 = XN+1 + PN+2 Te 

(1) 

The classical stochasticity parameter = ke. This map 
depends also on a further parameter — {T — e)/e. 

Clearly, we see that setting = 1 in Eq.([T|) recovers 
the Standard Map or DKP. To obtain the 2-DKP, we 
take ^ 10 — 100: typical experimental values used in 
are r, « 10 - 25 and K, = 0.1 to 0.5. The RP-DKP is 
obtained by taking the limit — » c», causing the impulse 
sinxjv+2 to become randomized. This can be achieved in 
practice by keeping k and e constant and taking T — > oo. 
In the usual T-scaled map this would yield an infinite 
stochasticity constant K — kT (but little insight). This 
RP-DKP limit has no mixed phase space behavior at all, 
but it retains the momentum trapping and its behavior is 
determined by the stochasticity parameter K,: . The map 
Eq.(IT]) has 27r periodicity in these re-scaled momentum 
units for integer r^. 

FiglD (a) compares Poincare surfaces of section (SOS) 
for the 2-DKP and RP-DKP. For the 2-DKP, a peri- 
odic structure of chaotic 'cells' separated by thin mixed- 
phase regions is apparent. These momentum 'trapping' 
regions appear at p ~ (2m -I- l)7r where to = 0, ±1, ±2.... 
For odd-integer multiples of tt there is near-cancellation 
of consecutive kicks. Fig[l] (b) shows the RP-DKP, 
for which — 10^. In fact, an indistiguishable SOS 
can be obtained by taking actual random phases, i.e. 
xn+2 = 27rCjv where < < 1 is a random number 
chosen every kick pair. We see that the trapping is there, 
but with no trace of islands etc. The RP-DKP never has 
any islands, regardless of how small K,: becomes. 

For the quantum equivalent we note that in a basis of 
plane waves, the one-period time-evolution operator for 



the 2-DKP, W^, has matrix elements 
where 



rrfree j-T2kick 



C^fm = J2 -^l-^ J^-^ exp {-iehj2) (2) 



where J„ is the integer Bessel functions of the first kind, 
?ie = he is the effective value of Planck's constant. Kn = 



2^ 1--^-> "-|- ! 
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FIG. 1: Comparison of the Poincare SOS for the 2-DKP and 
Randomised-Pair DKP. Stochasticity parameter = 0.4. 
(Left) 2-DKP: = 24. A periodic structure of chaotic 'cells' 
separated by thin mixed-phase regions (coinciding with the 
trapping regions) is apparent. (Right) RP-DKP: Tc = 10^ 
so impulses between kick pairs are uncorrelated. A similar 
plot is obtained by taking xn+2 to be a random number. The 
strength of trapping is largely unchanged, but all remnants of 
mixed-phase space structures are eliminated. Our key finding 
here is that the fractional scaling of the quantum localization 
lengths remains unchanged for the RP-DKP. 



an unimportant phase, it can be shown that [lo| : 
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2kick 
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COS [lh^/2] 



(3) 



aI—ttl 



exp i^^Te?ie/2) . To within 



Eq|3] can be derived most straightforwardly by evaluat- 
ing the sum in Eq[2] using Poisson summation. The key 
point is that Uf^'^^ is common to all the kicked sys- 
tems: only U(J^'^{t^) determines whether we investigate 
the standard DKP {t^ = 1), the mixed phase-space 2- 
DKP (1 < r, < 100) or the random RP-DKP [t, oo). 
Equivalently, we can also obtain the RP-DKP by taking 
Jjfree _ g^p ( — i27r^/) , i.e. using random phases for each 
angular momentum I. Taking = 10^ — 10^ gives the 
same behavior, provided is not a rational multiple of 
tt: transport can be strongly affected at such resonances. 

Figlllja) shows the asymptotic quantum momentum 
distributions, iV(p, t — > oo), obtained from an initial state 
il}{t = 0) '-^ 5{p — po) for pq = 0, in regimes analogous 
to FigllJ The dashed line corresponds to a momentum 
distribution obtained for the mixed phase-space 2-DKP 
while the solid line represents the equivalent result ob- 
tained using random phases in U^^'"^ . One sees that both 
the mixed phase-space 2-DKP and the random map (RP- 
DKP) distributions are very similar, with a distinctive 
'staircase' structure. 

Figllljb) shows the staircases obtained for — 0.4, 
% = 1/25,1/250 and Ti, = 1/2500 for the 2-DKP. We 
have rescaled the momentum distributions N{p) using 
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FIG. 2: (a) Shows that the quantum momentum distributions 
for the mixed phase-space double kicked atoms (2-DKP) are 
similar to those of the RP-DKP, for which the phases be- 
tween kick pairs are random numbers. The staircase struc- 
ture and 'heights' of the steps are comparable. Tfe — 0.4 and 
he = 1/500 in every case, (b) Shows the fractional scaling 
for the 2-DKP. = 0.4 and = 49. The JiJ'^ rescaling of 
the momentum distributions shows that the first step scales 
near-perfectly in the range from fi^ — 1/2500 to 1/25; i.e. the 
'local' localization length in the trapping region scales over 
this entire range. The global localization length (the enve- 
lope of the whole distribution) shows this fractional he scaling 
holds only for the smaller h^ = 1/2500 and he = 1/250 but 
not for he = 1/25. Vertical lines indicate width and centre 
{p = tt) of the first trapping region. The inset shows the same 
curves without the h^^'^ scaling. 
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FIG. 3: (a) Fractional scaling of the Ljoc (the localization 
length at the first step: note that Lioc ~ (2d)~^). It shows 
the fractional 2d ^ fi'J'^ scaling holds over the full he range. 
Straight lines indicate slope 0.65. The « 2/3 scaling holds 
even for Ke — 0.97 « Kc, the parameter for criticality in the 
Standard Map: the only difference is that here Te — > oo (RP- 
DKP) instead of = 1 (Standard Map), (b) As for (a) but 



shows the average of steps 2 to 4, showing that the L 
global scaling persists only for the smaller he values. 
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the 'local' scaling exponent estimated in Fig[3]for = 
0.4 which yielded a « 0.65 . We plot ln(7V(p));i"°■^^ 
shifted by an appropriate constant. The figure shows 
that the shape and magnitude of the distribution around 
the first trapping region re-scales perfectly with h^^ . Lioc 
is the slope on the first step. The height of the steps 
(2d in the notation of 10]) scales fractionally with h^, 
while the width w (momentum width between the dashed 
vertical lines) is independent of h^: in scaled momentum 
units w ^ 2tt/6. Thus we take Lioc ~ w/2d, i.e. for the 
scaling, Lioc ^ (2d)~^ 

The longer-ranged localization length L, characterizing 
the envelope of the full staircase, scales well in the smaller 
range < 1/250 but we note that the = 1/25 
staircase does not scale at long-range. 

In FigO we investigate the scaling exponent itself. 
For staircases where ?ie is sufficiently small, \n{N{p)) « 
const outside the 'steps' in the trapping regions. Then, 
only the step heights depend on . We can restrict our- 
selves to a study of the parameter 2c? since it has the 
same h^^ scaling as Lioc- The midpoint of the n-th step 
is at p^^\ — 2{n — l)7r/e for n 1; for n = 1, we take 

Pi/2 " '"'/(^c)- For the 7i-th step we measure 

2d(n) = (;n(7V(p(;^))) - (;n(7V(p(;+^)))) (4) 
where the average is taken over a small momentum in- 



terval around p\^2- ^^^^ procedure is easily automated, 
allowing a very fine-grid of ?ie values, so we can examine 
fiuctuations in behavior. 

FiglSlJa) compares the ?ie scaling of 2d{l) for the 2- 
DKP and the RP-DKP, while FiglStb) shows the aver- 
age of 2d for steps n~2 to 4. We take the behavior of the 
first step to indicate the local localization length and of 
steps 2 to 4 to indicate the properties of the long-ranged 
scaling. All the plots, whether 2-DKP or RP-DKP, show 
a slope a ~ 0.65 for small h^. For the upper ranges of 
h^, the slope (for L only) increases significantly. These 
results are quite consistent with the slightly higher esti- 
mate cr « 0.75 found in [13]: there, an average of the first 
few steps was obtained, including a few h^^ values in the 
range where the deviation from cr « 2/3 begins. 

We note that the trapping regions and the Standard 
Map golden cantori regions occupy similar regions of 
phase-space: the golden tori occur at winding number 
Po/27r « 0.618 and po/27r « (1 - 0.618). The trapping 
occurs for pq/2'k « 0.5 (e.g. momentum p « tt), ex- 
actly midway between these. In Fig.lU) we show the 
well-known plot of Geisel et al (Fig. 2(b) in Q) repro- 
duced in quantum chaos textbooks [llj. In that study, 
an initial state 'ip{t = 0) ~ 5{p — po) was evolved for crit- 
ical Kc ~ 0.9716 to study transport through the golden 
cantori. Notably, that study used precisely po = 3.2 tt: 
the initial state was centred in the trapping region. We 
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FIG. 4: Figure shows that the RP-DKP trapping regions 
coexist with the phase space regions which contain the two 
'golden ratio' tori of the DKP (Standard Map); nevertheless, 
destroying the cantori in the RP-DKP with random phases 
does not eliminate fractional scaling. (Bold solid lines) DKP 
at critical Kc ~ 0.97; the plot is similar to Fig. 2b of Geisel et 
al 1986 01 which yielded Lioc ~ K^^^'^ . Vertical lines indicate 
'golden' tori. Graph was obtained by setting = Kc = 0.97, 
= 1,K = 0.01. (Dotted line) the RP-DKP distribution is 
shown for a kick strength — 0.1, sufficiently small that the 
amplitude 'drops' are comparable in magnitude to the DKP 
(but without the golden cantori). 

compare a RP-DKP result also with po — 3.2 « tt but de- 
stroying all the KAM structures including cantori by tak- 
ing Te ^ oo but otherwise running exactly the same nu- 
merical code (We take the same h = 1/100 but K = 0.1 
for the RP-DKP to keep diffusion rates comparable to 
the DKP). Since e is our unit of time in both systems, 
we can drop the e subscripts on h and K. 

We now analyze the cause of the L ~ h^"^^^ scaling of 
the RP-DKP. We consider the trapping regions 5p ^ {p — 
Po) centered at po = IqTl = {2n + l)7r. At the centers of 
these regions the argument of the Bessel function in Eq[3] 
vanishes. Near these points it is of appreciable size up 
until when its argument x = ^ cos [111/2] is of the order 
of the index; i.e. when x \l — m\. If also |? — to| 3> 1, 
the asymptotic approximation 

t/fj;-^ ~ Ji.m{x) ^ ,-l!l^Mi-2'^'z), (5) 

with X = |Z — m| -t- |j_^|i/3 , holds. Sp is small, as we 
are near the center of a bottleneck, so this implies we are 
in the semiclassical, ?i — > 0, regime. Heuristically, the 
transition probability corresponding to a small change 
in momentum Ap is jC/;™*'^'^!^, which is proportional to 
|; _ to|-2/3 ^ (where Ap = \l - m|a). The effect 

of U^^'^'^ is to randomize the phases. The trapping re- 
gions act as bottlenecks for transport. On scales larger 
than Ap, considered fixed, diffusion takes place: p^ ~ Dt, 
where the local diffusion coefficient D is proportional to 
t^im*'^'^!^ and therefore to h^^"^. Assuming that the evo- 
lution operator is a band-matrix with constant width (as 



is the situation for the ordinary Kicked- Rotor), the lo- 
cal localization length satisfies 0, [l3| Liocai ~ ™" 
plying Liocai ^ fi ■ In our case, the width of the 
contributing region, as implied by the scaling of the ar- 
gument of the Airy function in Eq. (6), is of the order 
|/ — m\^/^ = ^h}/3 ■ This should multiply the expres- 
sions for D and Liocai- The resulting localization length 
therefore satisfies L ~ hr"^^^ . The condition for the va- 
lidity of the asymptotic properties of the Bessel function 
and for the momentum to be in a bottleneck region is 
1 > cos [in/2] > 

We recall that this structure of Uf^^^ is common to 
both the quantum DKP as well as RP-DKP; and FigH 
shows that the fractional scaling, is observed over similar 
phase-space regions. One may safely conclude from the 
analysis here that the behaviour of the generic RP-DKP 
follows only from the semiclassical dynamics in the bot- 
tleneck regions: it is the hallmark of the scaling of an Airy 
function rather than of the fractal scaling near golden 
cantori. However, for the DKP, the additional presence 
of large stable islands bordering the regions where 'h?/'^ is 
observed (which occur when « 1) make the quantum- 
classical interplay less transparent. In this case too. Airy 
functions and h^^^ factors occur naturally in the semi- 
classical quantization of torus states [l^]- A much more 
detailed study than the one undertaken here is required 
to conclusively establish whether the semiclassical dy- 
namics is also the dominant mechanism in the quantum 
DKP, but we suggest that cantori do not represent the 
only possible source for K^'^^'^ scaling behavior. 
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